SNSN-xxx-xx 
January 13, 2011 



Studies of exclusive charmless semileptonique B decays 
and extraction of \Vub\ at BaBar 



Martin Simare0 

Department of Physics 
Universite de Montreal, Montreal, CANADA 



We report on recent measurements of the branching fractions (BFs) for 
the decay channels S° 7r-£+i/, S+ 7r°£+i/, S+ r/£+z/, 5+ t/'^+z/, 
5° — )■ p~i~^h' and i?"*" — )• p^d'^y. We obtain very precise values of the total 
branching fractions for these decays, as well as partial branching fractions 
as a function of for the decay channels — 7r~£+z/, — )■ r](.'^v and 
B^ —7- p^l^u. In particular, we use the partial branching fractions of 
the 5° — )■ TT~i~^i> decay channel and form-factor calculations to extract 
several values of \Vub\- When we compared these values of \Vub\ to the one 
measured in inclusive semileptonic B decay, we find that two of them are 
consistent, within large theoretical uncertainties. 
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1 Introduction 



Semileptonic decays are best to measure \Vub\ P because they are much easier to un- 
derstand theoretically than hadronic decays and they are far more abundant than lep- 
tonic decays. In particular, the decay rate of B irii' is proportional to \ Vubf+{q'^)\'^, 
where f+{q^) is the theoretical form-factor calculation as a function of q^, the mo- 
mentum transferred squared. A precise value of \Vub\ from the exclusive B — )■ tiIp 
decay can be performed to test the QCD calculations and to constrain the description 
of weak interactions and CP violation in the Standard Model. 

We present two recent analyses in BABMi where \Vuh\ is obtained in the study of 
the exclusive B — )■ ixlv decay. In the tt — 77 analysis [2] , we study three decay modes: 
B^ —7- 7r^£+z/, B^ — 7- r}l^p and 5+ — t- rj'l^v. In the vr — p analysis [3], we study four 
decay modes: S° n-i+u, B+ n^i+u, 5° p'l+p and B+ p^l+p. These 
decays involve a 6 — )■ m transition via the coupling oi &W gauge boson. The value 
of \Vub\ can be extracted from the partial branching fraction, Ai3(g^), measured as a 
function of q^: 
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where r^o = 1.525 ± 0.009 ps |1] is the B^ lifetime. As can be seen in Eq. [H A(^(g^) 
depends on the f+{q^Y form factor that is provided by QCD calculations from light 
cone sum rules at low q^ values and lattice QCD at high q^ values. 



2 Experimental method 

The main differences between the n — r] and the tc — p analyses are summarized in 
Table m Both analyses use an untagged technique. This means that only one of the B 
mesons of the BB pair is reconstructed. Values of q"^ are determined using different 
methods which do lead to some variations in their values. The two-dimensional 
distribution of true versus reconstructed values of yields a detector response matrix 
which is used to unfold the measured distribution onto the true one. 

Backgrounds can be broadly grouped into three main categories: decays arising 
from b — )■ uiu transitions, decays in other BB events and decays in continuum events. 
Given the sufficient number of events in the tt — t] analysis for the 7T~i~^i> decay mode, 
the data samples can be subdivided in 12 bins of q^ for the signal. For the tt — p 
analysis and the other decay modes, a smaller number of events leads us to restrict 
the number of bins used in the fit. We use the AE-m-Es histograms, obtained from 
the Monte Carlo (MC) simulation, as two-dimensional probability density functions 
(PDFs), to extract the yields of the signal and backgrounds as a function of g^ in our 
fit to the data. 
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Analysis 

Luminosity on T(45') peak 

Number of BB pair events 

q'^ evaluation 

Cut strategy 

Cut selection 

Signal efficiency 

Background / signal 

5° ^ 7T-i+u yield 

Number of bins in vr mode 

Systematic uncertainties 



vr 



V 



422.6 fb~^ 
464 millions 

(.-Pb Pmeson) 

cuts, q^ dependent 
loose u cuts 
8% to 15% 

11.5 
11778 ±435 
12 

full gaussian 



TV — p 

349.0 fb-^ 
377 millions 

(Pi + Puf 

NN, q^ dependent 
tighter u cuts 
6% to 7% 
6.3 

10604 ± 376 
6 

±la 



Table 1: Comparison of various characteristics for the two recent analyses in BABML. 



In each analysis, the systematic uncertainties are estimated from the variations of 
the resulting partial BF values when the data are re-analyzed with different simula- 
tion parameters. In the tt — 77 analysis, for each parameter, we produce new PDFs by 
varying randomly the parameter value over a complete gaussian distribution whose 
standard deviation is given by the parameter uncertainty. One hundred such vari- 
ations are done for each parameter. The systematic uncertainty of a parameter is 
given by its RMS value of the resulting partial BF distribution from these one hun- 
dred variations. In the vr — p analysis, the systematic uncertainties are evaluated by 
±lo" variation for each parameter. 



3 Results 

The experimental AB{q^) distributions for B^ tt'I^u decays are displayed in 
Fig. [1] in the tt — 77 analysis, together with two parametrizations and three QCD 
calculations, and in Fig |2] in the tt — p analysis, combining the charged and neutral 
pion channels assuming isospin symmetry. From the BGL expansion extrapolated 
to q^ = 0, we obtain the value of |Ki6/+(0)| in both analysis as shown in Table [2l 
These values differs from both analysis since the experimental distributions do look 
different. However, the individual values of the partial branching fractions are indeed 
consistent with each other for the two analyses. The comparison between theory and 
experiment in their ranges of validity shows that all three QCD calculations are 
compatible with the data as we can see in Table [2l 
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Figure 1: Partial AS(g^) spectrum in 12 
bins of for — )■ ti^I^p decays in 
the n — rj analysis. The solid green and 
black curves show the result of the fit 
to the data of the BK [5J and BGL [6J 
parametrizations, respectively. The data 
are also compared to unquenched LQCD 
calculations (HPQCD [7], FNAL [H]) and 
an LCSR calculation [9]. 



Figure 2: Partial B — > tt^z/ spectrum 
in 6 bins of for B — t- n£iy decays in 
the 11 — p analysis. The red curve rep- 
resents the simultaneous fit to the data 
points and four theoretical points pro- 
duced with the FNAL LQCD calculation 
(magenta, closed triangles). 



Since, in the n — p analysis, the number of data points is limited to two above = 
16 GeV^, it was deemed desirable to undertake a simultaneous fit of theoretical and ex- 
perimental points to extract a value of \Vub\- We obtain the total BFs B{B^ — )■ 7r~£+z/) 
= {lA2±0.05stat±0.07syst) X 10-^ B{B+ ^ 7]£+u) = (0.36 ± 0.05,*,* ± 0.04,^,,*) x 
10"^ and B{B+ ri'i+u) = (0.24 ± 0.08,^* ± 0.03,j,,*) x 10"^ in the vr - r/ analysis, 
and B{B^ n-i+u) = (1.41 ± 0.05,^* ± 0.07,^,,*) x 10"^ and B{B° p'l+v) = 
(1.75 ± O.lbstat ± 0.27syst) X 10""^ in the tv — p analysis. Values of \Vub\ obtained in 
our two analyses are given in Table [2J They range from (3.0 — 3.8) x 10^'^. 



4 Summary 

It is estimated that there is less than 20% overlap in the selected event samples 
between the two analyses for the B^ — )■ -K'i'^iy decay channel. It is thus very satisfying 
to note that there is excellent agreement between the results of the two analyses. The 
values of the total BFs obtained in our work are the most precise total BFs to date. 
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Analysis n — r] vr — p 
HPQCD [7\ {q' > 16 GeV^ 3.24 ± 0.13 ± 0.161^;;^^ 3.21 ± 0.171^;^^ 
FNAL 18J {q^ > 16 GeV^) 3.14 ± 0.12 ± 0.16t|]:i 2.95 ± 0.31 
LCSR [9j (g2 < 12 GeV^) 3.70 ± 0.07 ± Omtoit 3.78 ± 0.13l[;i^ 
\Vubf +{0)1 (8.6 ± 0.3stat ± O.Ssyst) X 10"^ (10.8 ± 0.6) x 10"^ 



Table 2: Values of \Vub\ x 10 ^ derived from the form-factor calculations for the 
— i- n^i'^u decays. 

Our value of the total BF for — )■ rj'i^u, with a significance of 3. Oct, is an order 
of magnitude smaller than the CLEO result jlO]. The three values of \Vub\ are all 
acceptable according to the data. Two of them [71 |9] are consistent, within large 
theoretical uncertainties, with the value measured in inclusive semileptonic B decays: 
\Vub\ = (4.27 ±0.38) X 10-3 g]. 
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